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ABSTRACT 

We present 45^^-47'^ angular resolution maps at 50^' sampling of the ^^CO 
and ^^CO J=l-0 emission toward a 1.39° x 1.33° region in the W51 HII 
region complex. These data permit the spatial and kinematic separation of 
several spectral features observed along the line of sight to W51, and establish 
the presence of a massive (1.2 x 10^ ^o), large (A-^ x Ab = 83 pc x 114pc) 
giant molecular cloud (GMC), defined as the W51 GMC, centered at 
{i,b,V)c ^ (49.5°, -0.2°, 61 kms-^). A second massive (1.9 x 10^ M©), 
elongated (136 pc x 22 pc) molecular cloud is found at velocities of ^68kms~^ 
along the southern edge of the W51 GMC. Of the five radio continuum 
sources that classically define the W51 region, the brightest source at A 6 cm 
(G49. 5-0.4) is spatially and kinematically coincident with the W51 GMC and 
three (G48.9-0.3, G49. 1-0.4, and G49. 2-0.4) are associated with the 68kms-^ 
cloud. Published absorption line spectra indicate that the fifth prominent 
continuum source (G49.4-0.3) is located behind the W51 molecular cloud. The 
W51 GMC is among the upper 1 % of clouds in the Galactic disk by size and the 
upper 5-10% by mass. While the W51 GMC is larger and more massive than 
any nearby molecular cloud, the average H2 column density is not unusual given 
its size and the mean H2 volume density is comparable to that in nearby clouds. 
The W51 GMC is also similar to other clouds in that most of the molecular mass 
is contained in a diflFuse envelope that is not currently forming massive stars. 
We speculate that much of the massive star formation activity in this region has 
resulted from a collision between the 68kms~^ cloud and the W51 GMC. 
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Subject headings: ISM: clouds — ISM: general — ISM: individual (W51) — 
ISM: molecules 



1. Introduction 



The compact radio continuum sources comprising W51 ([Westerhout 19581) have long 



been recognized to constitute some of the most luminous star formation regions in the 
disk of the Galaxy ( [Harper fc Low 1971] ; [Hottman, l^rederick, fc Emery 1971| ). The high 



luminosity, the large number of inferred O type stars ( [Bieging 1975| ), and the location of 



these sources within a molecular cloud ( [Mufson fc Liszt 1979|) all suggest that the W51 



region represents the early formation stages of an OB association. Besides the intrinsic 
interest in the properties of W51, this region represents one of the closest analogs in the 
disk of the Milky Way to the more vigorous star forming sites found in other galaxies (e.g. 
30 Doradus). Since these latter regions are quite distant, W51 affords many advantages in 
investigating the detailed properties of luminous star forming sites and inferring how these 
regions may originate. 

One key to understanding the formation and evolution of any star forming region 
is establishing the properties of the molecular cloud out of which the stars form. While 
the molecular gas in the W51 region has been the subject of numerous studies, the 
interpretation of the results remain controversial. Scoville & Solomon (1973), primarily on 
the basis of small strip maps in CO (1-0), identified several molecular line components 
toward W51 and derived a minimum mass of 10^ Mq and a diameter ~ 20-30 pc for the 
molecular cloud that they associated with the most intense radio component at A 6 cm 
(G49. 5-0.4; [Mehringer 1994| ). They further suggested that this cloud might be physically 



related to the several thermal radio continuum sources that make up the W51 Hll-region 
complex ( IGoss fc Shaver 1970| ; [Wilson et al. 1970| ; [Georgelin fc Georgelin 197^ ) . Subsequent 



studies of the molecular gas toward W51 have confirmed the existence of a large molecular 
cloud ([Mufson fc Liszt 1979| ; [Nakamura et al. 1984| ; [Ohishi et al. 1984[ ; Pame et al. 1987] ; 



Solomon et al. 19"8^ ; [Scoville et al. 198^ ), although various models continue to be proposed 



for the relationship of the multiple spectral features seen in the molecular gas lines and 
their association with the different HII regions. 

The primary difficulty in establishing a definitive model of this region is the unique 
location of W51 in the Galaxy with respect to the sun. The W51 region has classically 
been associated with the tangent point of the Sagittarius spiral arm, which is oriented such 
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that the line of sight toward W51 intersects the spiral arm over several kiloparsecs of path 
length ( ^hane fc Bieger-Smith 196^ ; [Burton 1970|) . Much of the uncertainty surrounding 
the W51 region stems from establishing whether the numerous radio continuum sources 
and molecular clouds represent a single, large massive star forming region, or the chance 
projection of unrelated star forming sites viewed down a spiral arm. To better place 
the W51 region in context with respect to its location in the Galactic plane. Figure [l| 
displays the integrated ^^CO(l-O) emission in lOkms"^ velocity bins covering longitudes 
40°-55° from the Massachusetts-Stony Brook ^^CO Survey ( ganders et al. 198^ ). The W51 
region is distinguished by bright ^^CO emission extending over a l°xl° area centered on 
{i,b) - (49.5°, -0.2°) at velocities >55kms-^ A "3^D" view of the {i,b,V) ^^CO data 
cube covering the region surrounding W51 is shown in Figure |2|. The ^^CO isointensity 
contour surface in this figure clearly illustrates both the relatively large number of smaller 
molecular clouds with typical internal velocity dispersions of AV ~3-5 kms~^, and the 
large concentration of ^^CO emission extending over a ^20kms~^ interval in the W51 
region. Much of the ^^CO emission in this area has centroid velocities that exceed the 



maximum velocity permitted by pure circular rotation (i.e. 1/max ^54-57 km s ^; [Brand 



Blitz 1993| ). Such velocities have long been noted in 21 cm HI surveys at longitudes near 
i = 50°, and have been attributed to large-scale streaming motions of gas in a spiral density 
wave (e.g. [Shane fc Bieger-Smith 1966[ ; [Burton 1970| ). 

In principle the extent and properties of the molecular clouds located in the W51 region 
can be established by using the kinematic information in the molecular line data to isolate 
individual clouds. In practice, previous surveys have had either poor resolution or sparse 
sampling to make such an attempt feasible. Therefore, we have obtained full beam sampled 
maps of the W51 region in both CO (1-0) and CO (1-0) at subarcminute resolution in 
order to determine the relationship between the various molecular components. These maps 
permit us to disentangle the blends of unrelated clouds along the line of sight and to obtain 
more accurate mass estimates of the molecular gas. These data can also be compared with 
similar maps of more nearby clouds that have recently been obtained by us and others. 

The outline of this paper is as follows. In Section ||, the observing procedures are 
described and channels maps of the ^^CO and ^^CO emission are presented. Analysis of 
the diflPerent spectral features observed in these maps and a more thorough discussion of 
the main features associated with the compact radio continuum sources in W51 is given in 
Section |^. In Section ||, we discuss the current massive star formation events in the region 
with respect to the various molecular components and comment on the evolution of the 
W51 Giant Molecular Cloud (CMC). Our conclusions are summarized in Section |. 
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2. Observations 

2.1. Observing Procedure 

A 1.39° X 1.33° region (100 x 96 pixels) toward the W51 region was mapped in 
i2CO(l-0) (115.271203 GHz) and ^^00(1-0) (110.201370 GHz) in April 1995 using the 
QUARRY receiver array ([Erickson et al. 1992|) on the FCRAO 14 m telescope. The fuU 
width at half maximum beam size of the 14 meter antenna at these frequencies is 45^' 
and 47^' at 110 and 115 GHz respectively. The data were taken in position switching 
mode and calibrated with the standard chopper wheel method of observing an ambient 
temperature load and sky emission. The backends for each pixel of the array consisted 
of an autocorrelator spectrometer set to span the velocity range from 0-100 km s~^ at 
78 kHz sampling (0.20 km s~^ @ 115 GHz) and 94 kHz resolution (0.24 km s~^). During 
data reduction the spectra were smoothed to a velocity resolution of 0.5kms~^. 

Previous FCRAO measurements indicate that the spillover and scatter efficiency (^s) 
of the telescope and radome is ^0.7 at the observed frequencies. The observed antenna 
temperatures corrected by ryfss are presented as ( [Kutner fc Ulich 1981|) . A further 



correction, the source coupling efficiency (rye), accounts for the coupling of the beam to the 
source. For a uniform source that fills only the main beam of the 14 m telescope, r]c is ~ 0.7 
(i.e. 30% of the power is scattered on angular sizes much greater than the FWHM beam 
size), while for sources with uniform intensity over a diameter of 30', rjc is 1.0. In practice, 
the observed structures in the ^^CO and ^^CO maps span a range of sizes and shapes, and 
applying a single coupling efficiency for the entire map is incorrect. For simplicity, we 
present and analyze the data in the temperature scale. The typical rms noise in the 
^^CO and ^^CO maps in 0.5kms~^ channels is AT^ ^0.7 K and 0.6 K respectively. 



2.2. The Data 

Images of the integrated ^^CO and ^^CO intensity (/T^ dv) in 2kms~^ wide intervals 
are presented in Figures ^ and |^ respectively in the velocity range from 40 to 70kms~^. 
The values printed in the upper left corner of each figure panel denote the centroid velocity 
of the particular interval. Extended ^^CO and ^^CO emission was detected between Okms~^ 
and 25kms~^, but these data are not presented here. This low velocity emission most 
likely originates from local molecular clouds and is not related to the W51 region of interest 
here. Little emission was observed between 25 and 35kms~^ and at velocities in excess of 
75kms~-^ (see §|). Similar velocity structure is also observed in the 21 cm HI emission 
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lines ([Burton 1970| ). The following section analyzes the velocity structure in the molecular 



line maps and identifies individual molecular clouds. 

3. Analysis 

3.1. Velocity Structure 

The ^^CO and ^^CO emission toward the W51 region contains a number of discrete 
velocity components that overlap in projection both spatially and kinematically. To identify 
and isolate the emission from these velocity components, multiple gaussians convolved with 
the spectrometer channel widths were fitted to each spectrum in an automated manner. The 
free parameters for each gaussian were the amplitude of the spectral line, the mean velocity, 
and the line width. The number of gaussians fitted to each spectrum was determined by 
searching for contiguous channels that contain an integrated intensity with a signal to noise 
ratio > 3. Channels containing a local antenna temperature maxima (denoted here as 
channels q, i = 1, n) in each such section were then identified. A local maximum at channel 
Ci was deemed a "significant" peak if the antenna temperature in any channel between c^ 
and the neighboring local maximum at channels q_i and q+i decreased by more than 2 arms 
from the antenna temperature at channel q. The number of significant peaks corresponded 
to the number of gaussians fitted to that section. Each spectrum was smoothed to a velocity 
resolution of 1.5kms~^ prior to identifying the peaks, although the fits were performed on 
the 0.5kms~^ resolution data. Spectra with large residuals with respect to the gaussian 
fits were visually inspected and additional gaussians were added as appropriate. The ^^CO 
data were easily decomposed into gaussians in this manner, but it often became difficult 
to reliably identify the velocity features in the heavily blended ^^CO lines. Also, toward 
the compact HII regions, some of the structure in the ^^CO and ^^CO line profiles can 
be attributed to the absorption of radiation from hot molecular gas by colder foreground 
material (see [Mutson fc Liszt 1979|) . Away from these compact regions, absorption effects 



are not as significant, and over most of the cloud, the peaks in the spectral lines should 
accurately represent the velocity structure along the line of sight. 

The results from the gaussian decomposition of the line profiles are synthesized in 
Figure]^. The upper panel shows histograms of the mean velocities for the ^^CO (thick lines) 
and ^^CO (thin lines) gaussians, and the lower panel shows the total integrated intensity 
in the gaussians as a function of the mean velocity. Most of the emission is confined to 
velocity intervals of 0-25 km s~^ and 35-75 km s~^. We identify the 0-25 km s~^ emission as 
originating from nearby molecular material and the 35^75 km s~^ emission with molecular 
gas in the Sagittarius spiral arm. Several velocity components occur repeatedly in both the 
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■^^CO and ^^CO gaussian fits as signified by the histogram peaks shown in the top panel in 
Figure |5|. In particular, velocity components at 7, 15-25, 44, 49, 53, 60, 63, and 68kms~^ 
are readily apparent. In terms of the ^^CO and ^^CO integrated intensity, the two major 
features are the 60kms~^ and 63kms~^ components. A 58kms~^ component is indicated 
as well since that is the centroid velocity of the molecular line emission toward the brightest 
radio continuum source in the W51 region QMufson fc Liszt 1979| ). Note, however, that the 
58kms~^ component is not a prominent feature as judged from Figure ^. The following 
discussion briefiy highlights the morphology of the individual velocity components. 



3.1.1. 7 and 15^25 km s ^ components 

The 7kms~^ velocity component contains weak, narrow lines over nearly the entire 
mapped region and is undoubtedly a nearby molecular cloud. The 15-25 km s~^ interval 
appears to contain a few distinct velocity features (see Fig. |^), but it is unclear whether or 
not these components are physically related. 



3.1.2. Ji.Ji.kms ^ component 

The 44kms~^ cloud is elongated parallel to the Galactic plane at 6 = 0°, although this 
cloud may form part of a larger structure that extends to lower Galactic latitudes. The 
molecular gas at these lower latitudes occurs at velocities of ^40kms~^, which is outside 
the velocity range assigned to this feature. 



3.1.3. 49 and 53 km s ^ components 

The ^^CO emission from the 49 and 53kms~^ components is more fragmented than 
the other features mentioned so far. These fragments may represent either individual 
clouds or the remnants of a once larger cloud in the Galactic plane. The 53 km s~^ cloud 
is distinguished by bright -^^CO emission near (-^, h) ~ (49.4°, —0.3°) that is associated with 
the compact HII region G49.4-0.3 (see §^T|). 
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3.1.4' 58, 60, and 63km s ^ components: The W51 Giant Molecular Cloud 

The 63kms~^ component extends for nearly a degree in length and is found mainly in 
the central and eastern part of the mapped region. This is best observed in the 66kms~^ 
panel shown in Figure which represents the line wing emission of this velocity component 
(as well as emission from the GSkms"-*^ cloud discussed below). The 60kms~-^ component 
consists predominantly of a diffuse patch of emission that extends into a filament to the 
east, and a second filament to the south. The 60kms~^ and 63kms~^ velocity components, 
along with the 68kms~^ cloud discussed below, likely correspond to the "high velocity 
stream" of 21 cm HI emission ([Burton 1970|) that has been attributed to the streaming 



motions of gas down the Sagittarius spiral arm. 

Careful inspection of the channel maps indicates that the spatial distribution of the 
60kms~^ and 63kms~^ components generally do not overlap. For example, the western 
edge of the 63kms~^ component closely matches the eastern edge of the 60kms~^ emission. 
This is best seen in Figure ^ and the 60kms~^ and 66kms~^ velocity panels in Figure |. 
Further, the extended emission from the eastern portion of the 63kms~^ component is just 
above the filamentary extension of the 60kms~^ component. Such interfaces are unlikely 
to occur by chance from two unrelated clouds along the line of sight, and suggest that the 
60 and 63kms~^ components represent kinematic structure within a single molecular cloud. 
Velocity differences of this magnitude are commonly observed in nearby molecular clouds 
(e.g. IBally et al. T987| ). 

The 58kms~^ component is dominated by bright, compact molecular emission and 
does not contain the diffuse extended emission that characterizes the 60 and 63kms~^ 
features. Inspection of the channel maps suggests that this velocity component also refiects 
the interval kinematic structure within a single cloud encompassing the 60kms~^ and 
63kms~^ clouds. For example, the emission from the filament protruding to the southern 
portion of the mapped region contains primarily a centroid velocity of 58kms~^ closest to 
bright compact ^^CO and ^^CO emission. Further away from this bright, compact emission 
region, the velocity of the filament changes to ^60kms~^. Similar velocity patterns are 
observed in emission features along the eastern and western edges of the map. These results 
suggest that the emission constituting the 58, 60, and 63kms~^ components represent the 
internal velocity structure within a single molecular cloud. Koo (1997) reached similar 
conclusions concerning the atomic hydrogen clouds at these velocities based upon HI 
absorption observations toward the radio continuum sources. Since the bright ^^CO 
emission associated with the 58 km s~^ and 60 km s~^ components are coincident with the 
brightest radio continuum source in the W51 region (G49.5-0.5; see §^.1|), we henceforth 
refer to the 58-60-63 km s~^ components as the W51 molecular cloud. 
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3.1.5. 68 km s ^ component 

The 68kms~^ component extends for ^ 1° east-west across the image and also hkely 
constitutes part of the "high velocity stream" identified in HI surveys ( [Burton 1970|) . 



Interestingly, the 68 km s~^ filament is located at the southern edge of the diffuse emission 
associated with the W51 molecular cloud at velocities ^63kms~^ (see Figs. | and ||). 
Again, such a clear truncation of the W51 molecular cloud at the location of the 68kms~^ 
cloud is unlikely to occur from two random clouds along the line of sight, and strongly 
suggests that these two clouds are physically related objects at a common distance. 
Nonetheless, the elongated appearance of the 68kms~^ cloud is in stark contrast to the 
roughly circular shape of the W51 cloud, indicating that these two objects are best treated 
as individual structures rather than a single molecular cloud. 



3.2. Properties of the W51 and the 68kms ^ Molecular Clouds 

The physical properties of the molecular clouds in the W51 region can be determined 
from the gaussian decomposition of the line profiles. We single out these two clouds since 
they contain four of the five bright HII regions found in radio continuum surveys (see 
^|0|). In deriving the properties, the distance to the W51 cloud is assumed to be 7.0 ± 
1.5 kpc as determined from proper motion studies of the W51MAIN H2O maser in the 
G49. 5-0.4 dense core (|Genzel et al. 19811) . The 68kms~^ cloud was assumed to have the 



same distance based on its apparent association with the W51 cloud as discussed above. 
The other clouds are not included in this analysis since their distances are unknown. 

The properties of the W51 and 68kms~^ clouds are summarized in Table [1|. The cloud 
size represents a visual estimate of the extent of the detectable ^^CO emission along the 
major (^max) and minor (^min) axis of the cloud. The cloud line width, A]/fwhm, is the full 
width at half maximum of the sum of the gaussian fits comprising the respective clouds. 
Two estimates of the cloud mass are provided in Table [l|. The virial mass, Mvir, was 
calculated using the expression Mvir = 209 R Al/^ M©, where /SV is the full width at half 
maximum line width in kilometers per second and R is the cloud radius (= \/^max^min/2.0) 
in parsecs at the zero intensity level. While the cloud size is actually measured at a finite 
^^CO intensity level, we find it unlikely that the clouds are appreciably larger at lower 
intensities, and no correction to the observed cloud size was applied. Note that the above 
expression for the virial mass is appropriate for a uniform density, spherical cloud. The 
equivalent expression for a and density cloud would decrease the constant factor in 
the virial mass expression to 188 and 125 respectively. 
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A second mass estimate can in principal be obtained from the ^^CO and ^^CO data 
using the LTE analysis ([Dickman 1978|) . However, in blended regions, it is often difficult 



to associate ^^CO gaussian fits with analogous ^^CO features. Therefore, the H2 column 
densities were estimated by applying a constant conversion factor to the ^^CO integrated 
intensities (pickman 1975| ; [Solomon et al. 1987| ; Strong et al. 1988| ). To ensure that the 



Galactic conversion factor is indeed valid for the W51 region, the conversion factor was 
estimated from lines of sight with unblended ^^CO and ^^CO lines in the 56-71 km s~^ 
velocity interval that defines the W51 and 68kms~^ clouds. The H2 column densities 
for these lines of sight were estimated using the procedure outlined by Dickman (1978) 
and assuming a ^^C0/H2 abundance of 1.5 x 10~^. A histogram of the ratio of the H2 
column densities to ^^CO integrated intensities is strongly peaked with a mean value 
of 1.7 X 10^° cm~^ (Kkms~^)~^, and is similar to the conversion factor that has been 
derived for the Galaxy (between ^ 2-3 x 10^^ cm" -2 (Kkms-i)-i; [Strong et al. 1988| and 



references therein). The masses (Mxco) computed from the ^^CO integrated intensities 
were calculated by adopting the conversion factor derived from the W51 data, and include 
a multiplicative factor of 1.36 to incorporate the mass contribution from heavier elements. 

The results summarized in Table |l| indicate that the W51 molecular cloud has a mean 
diameter of ^97 pc and a mas^ of ~ 10^ ^0- The similarity in the two mass estimates 
indicate the gravitational potential energy is approximately equal to the kinetic energy, and 
that self-gravity must play a critical role in the evolution of the W51 molecular cloud. The 
68kms~^ cloud is 136 pc in length but only 22 pc wide over most of the minor axis. Such 
an obvious departure from spherical symmetry renders the virial mass estimate suspect, and 
so the mass of the 68kms~^ cloud was estimated only using the ^^CO conversion factor. 
The mass obtained, ^ 10^ M©, is an order of magnitude less than that of the W51 cloud. 

A comparison of the W51 cloud properties with the size and mass spectrum of molecular 
clouds in the Galaxy (e.g. [Sanders, Scoville, fc Solomon 1985| ; [Solomon et al. 1987[ ) shows 



that the W51 molecular cloud is one the largest GMCs in the Galactic disk. Among the 
^ 5000 molecular clouds with diameters in excess of ^ 22 pc (and corresponding masses 
- 10^ Mq; i.e. GMCs), the W51 GMC is in the upper 1 % of the clouds by size and the 
upper 5-10% by mass. For the 68kms~^ cloud, although its mass is typical of a relatively 
low mass GMC, a distinguishing feature is its shape. The ratio of the major to minor axis 
of the 68kms~-^ cloud is ^6. In the cloud catalog by Solomon et al. (1987), 85% of the 



^The virial mass derived here agrees well with that estimated by Dame et al. (1987), but is a factor of 
4 less than found by Solomon et al. (1987). Solomon et al. (1987) derived a substantially larger velocity 
dispersion, and hence computed a larger virial mass, most likely by including the emission from the 68 km s~^ 
cloud and several other velocity components that we did not associate with the W51 molecular cloud. 
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clouds have an aspect ratio less than 2, and only one object has an aspect ratio greater 
than observed for the 68kms~^ cloud. Clearly an elongated shape over such a large length 
scale is unusual in the molecular interstellar material and may indicate that the 68kms~^ 
cloud is a transient structure originating from a relatively recent dynamical event. 



4. Discussion 

4.1. Massive Star Formation in the W51 Region 

Now that the individual molecular clouds in the W51 region have been identified, 
their relationship to the massive star forming sites can be explored. Three images of the 
W51 region are shown in Figure ||: a map of the integrated ^^CO(l-O) intensity generated 
from gaussian fits with mean velocities between 56 and 71kms~^ (i.e. the W51 and 
GSkms"-*^ clouds), an image of the 60/xm emission from the IRAS Sky Survey Atlas, and 
a A21 cm radio continuum image ( [Koo fc Moon 19971) . The bulk of the 60/xm emission is 



elongated parallel to, but slightly below, the galactic plane, and is coincident with bright 
radio continuum emission. Of the sources labeled in the A21 cm continuum map, W51C 
has predominantly a non-thermal continuum spectrum and is thought to be a supernova 
remnant (|Subrahmanyan fc Goos 1995|) , while the other sources have thermal spectra 
and are compact HII regions. The radio continuum sources G49.4-0.3 and G49. 5-0.4 are 
classically referred to as W51A (|Kundu fc Velusamy 1967| ), with the G49. 5-0.4 region 



containing the infrared source W51IRS1 ([Wynn- Williams, Becklin, fc Neugebauer 1974| ) 



and the H2O masers W51 North, W51 South, and W51MAIN ( |Genzel fc Downes 19771 ). 
Sources G48.9-0.3, G49. 1-0.4, and G49.2-0.4 are collectively known as W51B ([Kundu| 
fc Velusamy 19671) . The observed radio continuum fiuxes (|Koo 19971) and far-infrared 



luminosities (^10^ L©; [Rengarajan et al. 1984| ; [Harvey et al. 198^ ) imply the presence 



of one or more O stars in each of these regions. 

Figure |^ shows that many of the radio continuum sources have corresponding peaks in 
the molecular line and far-infrared images. In addition to the spatial coincidences, these 
radio continuum sources have recombination line velocities ([Wilson et al. 19701) similar to 



the velocity components identified from the ^^CO observations. The G49. 5-0.4 HII region 
has a recombination line velocity of 59kms~^ and is spatially coincident with the strong 
^^CO emission associated with the W51 molecular cloud. The recombination line velocities 
toward G48.9-0.3, G49. 1-0.4, and G49. 2-0.4 are 66, 72, and 66kms-^ respectively and 
are located along the 68kms~^ molecular cloud. Finally, the G49.4-0.3 HII region has a 
recombination line velocity of 53kms~^ and is coincident with the bright ^^CO and ^^CO 
emission from the 53kms~^ molecular cloud. (Note that the molecular gas associated with 
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this HII region is outside the velocity range that defines the W51 GMC.) HI (|Koo 1997 



and H2CO ( |Arnal fc Goss 1985| ) spectra toward this source exhibit absorption features 



at velocities of 63-65 km s~^, indicating that G49.4-0.3 must be located behind the W51 
GMC. However, these observations do not indicate whether this source is situated just 
beyond W51 and hence is physically related to the HII region complex, or if it is a distant, 
unrelated background massive star forming site. The molecular line maps presented here 
provide no compelling reason for (or against) such an association. 

To search for any additional massive star forming sites in the W51 molecular cloud, the 
IRAS Point Source catalog was examined for objects that have at least two "high" quality 
detections among the four IRAS bands and a rising spectral energy distribution toward 
longer wavelengths. These criteria were designed to select a reliable sample of objects with 
far-infrared colors characteristic of embedded star forming regions ( | Walker et al. 1055| ). 
Of the 40 IRAS point sources within the mapped region that meet these criteria, the 
seven brightest objects at 25/xm are located along the interface between the W51 and the 
68kms~^ clouds. Visual inspection of IRAS 60/xm image in Figure || confirms that the 
brightest sources are located along this ridge. Assuming that the other point sources are 
located at the distance of the W51 molecular cloud (although many of them are almost 
certainly foreground objects), point sources found away from this interface have far-infrared 
luminosities in the four IRAS bands less than 40,000 Lq and inferred spectral types later 
than ZAMS BO. 5 ( |Vacca, Garmany, fc Shull 1996|) . Thus embedded O type stars in the 



W51 region currently are confined to the southern extreme of the molecular cloud. 



4.2. Comparison to Other Molecular Clouds 

The extreme star formation characteristics of the W51 region raises the question as 
to whether the star formation activity stems from unusual global properties in the W51 
molecular cloud, or from unusual conditions found local to the massive star forming sites. 
These possibilities can be explored by comparing the W51 cloud properties with other star 
forming regions that have also been extensively mapped in ^^CO and ^^CO. In particular, we 
shall compare the W51 and 68kms~^ clouds with the molecular clouds associated with the 
HII regions Sh 140, Sh 155, Sh 235, Sh 247, Sh 252, and Sh 255 which have been mapped 
with the same receiver and telescope used for the W51 observations ( [Heyer, Carpenter, fc| 



Ladd 1996| ; [Carpenter, Snell, fc Schloerb 1995|) . The embedded high mass stellar content in 



this comparison sample is generally limited to a single early^B/late 0-type star and is in 
stark contrast to the cluster of O stars forming in the W51 molecular cloud. 
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The clouds in the comparison sample have diameters of ^ 20-55 pc and masses of 
^ 2 X 10^ to 7 X 10^ Mq. Thus the W51 molecular cloud is ^2.5 times larger and > 10 times 
more massive than these objects. Note that the molecular clouds associated with Sh 247, 
Sh 252, and Sh 255 are distinct regions within the Gem OBI cloud complex, which has a 
total mass of 3 x 10^ Mq and a diameter of ~ 150 pc ([Carpenter, Snell, fc Schloerb 1995| ). 



While the spatial size of the Gem OBI complex is larger than the W51 cloud, the W51 cloud 
exhibits a continuous structure ^ 100 pc in size as opposed to the fragmentary appearance 
of the Gem OBI complex. The 68kms~^ cloud is also more massive than the objects in the 
comparison sample and is significantly more elongated than any of the clouds considered 
here. Thus the W51 and 68kms~^ clouds are at the extreme in terms of cloud masses 
and sizes compared to objects in this sample. With the possible exception of the 68kms~^ 
cloud, these clouds are similar though in that they appear to be gravitationally bound. 

The size and mass of a cloud are not necessarily the key parameters that control the 
star formation activity within the cloud. Intuitively, one might expect that the amount of 
matter above a critical density to be the critical variable for otherwise similar clouds. Thus 
if the massive star formation activity in W51 has resulted from the large scale collapse of 
the cloud, one would expect the volume and column density densities to be larger than 
found in a typical cloud. As an indirect measure of the H2 column densities. Figure ^ 
shows histograms of the observed ^^CO integrated intensities for each of the clouds in our 
sample. The lowest integrated intensity shown for any cloud is 3 Kkms~^ since that is 
approximately the highest 3(7 detection limit among the various ^^CO surveys. Comparison 
of the ^^CO and ^^CO intensities suggests that the ^^CO emission is optically thin if the 
^^CO and ^^CO excitation temperatures are equal as assumed in the LTE analysis (see 
Dickman 1978|) . Therefore, the distribution of ^^CO integrated intensities should accurately 
trace the H2 column density distributions as along as the ^^CO abundance is roughly 
constant within a cloud. The ^^CO integrated intensities in Figure |^ can be converted 
to H2 column densities for an assumed ^^C0/H2 abundance of 1.5 x 10~^ ( [Bachiller fc 
Cernicharo 1986| ) with the formula 



A^(H2) = 3.05 X 10^^ Texe^-^^/^^^ j T^C^CO)dv cir-^. 

For an excitation temperature of Tgx = 10 K, the 3 Kkms~^ integrated intensity limit 
imposed for Figure |^ corresponds to an H2 column density of 1.5 x 10^^ cm~^, or an visual 
extinction of ^1.5"^ ( [Bohlin, Savage, fc Drake 1978|) . Variations in the ^^CO abundance 
and excitation conditions will obviously eflFect the absolute conversion from ^^CO integrated 
intensities to H2 column densities, and the comparisons here are intended to search for large 
differences (factors of several or more) in the typical column densities in these clouds. 



-13- 



Figure |^ shows that the distributions of ^^CO integrated intensities peak near the 
detection hmit of 3 Kkms~^ for each cloud, with a long tail toward the higher integrated 
intensities. The tail of these distributions correspond to the high column density regions 
and are often associated with star forming sites. The mean ^^CO integrated intensity among 
the clouds varies between 4.9 and 9.7 Kkms~^, with the W51 cloud containing the third 
highest mean intensity and the 68kms~^ cloud the highest. The high values found for the 
68kms~^ cloud may be a result of an unusual viewing angle, as either this cloud is a sheet 
of gas observed edge on or a long, narrow filament. These results imply that most of the 
mass in each cloud is contained in lines of sight with column densities corresponding to less 
than a few magnitudes of visual extinction. Thus the W51 GMC is similar to other clouds 
in that the diffuse envelope contains more mass than the high column density cores. 

While the W51 cloud contains a higher column density on average than the other 
clouds, this can attributed to the fact that it is more than twice as large as some clouds in 
the sample. Indeed, assuming that the 1.2 x 10^ Mq W51 cloud is distributed in a sphere 
of diameter of 97 pc (see Table |l|), the average H2 volume density is 40 cm~^, comparable 
to the volume density inferred in nearby molecular clouds ([Blitz 1991] ; [Carpenter, Snell,"^ 



Schloerb 1995[) . This suggests that the entire W51 molecular is probably not in an advanced 



stage of collapse, and that the intense star formation activity in W51 likely results from 
forces acting on a localized region. In retrospect, this is perhaps not surprising given that 
the massive star forming regions in W51 are located at the edge of the cloud and not in the 
center as expected if, for example, the entire cloud was systematically collapsing. 

Contrary to the global properties, the gas properties local to the W51 massive star 
forming regions do appear to be unusual compared to the molecular clouds in the solar 
neighborhood. Submillimeter continuum observations have shown that the core containing 
the G49. 5-0.4 HII region contains ^10^ Mq of gas with a mean H2 volume density of 
^5 X 10^ cm~^ over a 3 pc radius region ( [Sievers et al. 19"9TD . By contrast, submillimeter 
observations indicate that the most massive cores in nearby molecular clouds typically 
have masses ^2-3 orders of magnitude less than that of the G49. 5-0.4 core (e.g. jOldham 



et al. 19941 ; [Mezger et al. 1992[ ; [Schloerb, SneU, fc Schwartz 19871 ; [Jafte et al. r984[) . While 



the G49. 5-0.4 core does not necessarily contain higher gas densities, it does contain more 
gas at the densities needed to form stars. 



4.3. Evolution of the W51 Molecular Cloud 



The above discussion suggests that the key to understanding the massive star formation 
activity in the W51 complex is determining the forces that acted on a localized region 
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within the cloud. The molecular line maps presented here allow us to speculate on what 
these forces may be. As shown in Figures ^ and ^, the diffuse emission from the W51 GMC 
truncates at the location of the 68kms~^ cloud for velocities ^63kms~^. This morphology 
was used in §^ to argue that the W51 and 68kms~^ clouds are at a common distance since 
such an interface is unlikely to result from the chance superposition of unrelated molecular 
clouds. One way such an interface could form is if the 68kms~^ cloud has collided into the 
W51 GMC (see also |Arnal fc Goss 1985| ; [Pankonin, Payne, fc Terzian 1979| ). Molecular 
gas does extend below the 68kms~^ cloud at lower velocities, however. In this picture, 
given the three dimensional structure of the W51 GMC, this material has not crossed the 
path of the 68kms~^ cloud. Both H2CO and HI spectra toward the G49. 5-0.4 HII region 
( |Koo 199"7| ; |Arnal fc Goss 1985| ; [Mufson fc Liszt 1979|) contain absorption lines at velocities 
>65kms~^, and relative to the line of sight, the 68kms~^ cloud must be located in front 
of the W51 GMC. Since the velocity difference between the two clouds (at least 5 kms~^) 
is larger than the sound speed, a shock front will form that will compress the molecular gas 
and possibly induce star formation QElme^reen fc Elme^reen" 19781) . 

The qualitative model of two colliding clouds accounts for several properties of the 
W51 region. First, one would expect that star formation should occur preferentially along 
the interface region. The IRAS image in Figure ^ shows that this is indeed the case for 
massive stars, as the northern half of the W51 GMC appears devoid of embedded O type 
stars despite containing most of the cloud mass. The cloud collision model would also 
suggest that the massive star forming regions along the collision interface should have a 
common age. While the ages of these stars are not known, the lifetime of the mid-0 type 
stars found in the W51 region ([Mehringer 1994| ; see also [Koo 1997|) places an upper limit to 
the stellar ages of ^ 5 Myr (|Meynet et al. 199^) . The actual ages may be considerably less 
since these O stars are still in the compact HII region phase, which has a lifetime of less 
than 1 Myr (|Churchwell 1990| ; [Uomeron fc Torra 1996|) . Further support for a cloud-cloud 
collision model comes from considering the projected distance (73 pc) between the two 
furthest separated star forming regions along the W51/68kms~^ cloud interface. The 
sound travel time across this distance is nearly two orders of magnitude larger than the 
O star lifetime. Therefore, these star forming regions must have been created by either a 
single event operating along the entire southern edge of the cloud or from up to 4 separate, 
but nearly simultaneous, events. Given the scarcity of embedded O stars in the Galaxy, the 
cloud-collision model would provide a natural explanation for simultaneous star formation 
along the ridge. Ultimately, the suggested collision between the W51 and 68kms~^ clouds 
may be related to a spiral density wave. The anomalous gas velocities in the W51 region 
have long been attributed to streaming motions in the Sagittarius spiral arm (|Shane fc| 
Bieger-Smith 196^ ; [Burton 1970|) . The associated spiral density wave may have indirectly 
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led to the massive star formation activity in the W51 region by enhancing the number 
density of clouds and increasing the probabihty of a cloud-cloud collision. Further, a spiral 
wave shock, if present (see [Lubow, Balbus, fc Cowie 1986| ), will compress and flatten any 



clouds ( IWoodward 1977| ). Indeed, such a shock could account for the highly elongated 



shape of the 68kms ^ cloud. The W51 GMC, however, remains roughly circular in shape, 
and globally its evolution is likely still dominated by self gravity. 

Finally, we briefly consider the implications of these results for other star forming 
regions. While the W51 star forming region exceeds all nearby embedded star forming sites 
in terms of the number of O stars, bolometric luminosity, and dense core mass, W51 itself 
is dwarfed by some star forming regions in nearby galaxies. Most notably, the 30 Doradus 
region in the Large Magellanic cloud contains an order magnitude more O stars than W51 
( |Vacca et al. 191?5D . At larger distances, many interacting galaxies contain even yet more 



vigorous massive star forming regions that may be dense enough to represent young globular 
clusters ( [Whitmore et al. 1993|) . Most of the molecular gas appears to have been dispersed 
in these clusters already, and in any event, the distances to these systems precludes any 
detailed studies of the natal clouds. The most signiflcant piece of information aflForded by 
the W51 molecular maps is that despite containing one of the most massive dense cores 
known in our Galaxy, most of the mass in the W51 molecular cloud is not currently forming 
massive stars. Thus it is easy to imagine that the increase in the number of cloud-cloud 
collisions that presumably results in interacting galaxies may lead to a larger number of star 
forming regions throughout a single cloud or more intense star formation activity within 
a small region. In fact, the mass within the W51 molecular cloud is comparable to that 
in globular clusters, and it may not require the conglomeration of many giant molecular 
clouds to form these stellar systems, but the large scale collapse of a single GMC. 



5. Summary 

We have mapped a 1.39° x 1.33° region toward the W51 HII region complex at 45^^-47^' 
resolution and 50^' sampling in the J=l-0 transitions of ^^CO and ^^CO. From these data 
we have identifled the major molecular clouds and have associated these clouds with the 
massive embedded star forming sites in the W51 region. We flnd that: 

(1) The two most prominent clouds in the W51 region are the 58-60-63 km s~^ cloud 
(defined as the W51 GMC) and the 68kms-^ cloud. The W51 GMC is associated with the 
brightest A 6 cm continuum source in W51 (G49.5-0.4), and the 68kms~^ cloud contains 
the HII regions G48.9-0.3, G49. 1-0.4, and G49. 2-0.4. Published absorption line spectra 
( |Koo 19971 ; |Arnal fc Goss 1985|) indicate that the fifth bright H II region in the area. 
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G49.4-0.3, must be located behind the W51 GMC, but it remains unclear whether or not it 
is physically associated with the other star forming regions. 

(2) The mass of the W51 GMC and the GSkms"^ clouds are -1.2x10^ and 
1.9 X 10^ Mq respectively. The W51 molecular cloud is roughly circular in shape with 
a mean diameter of ^ 97 pc and appears to be gravitationally bound. Compared to the 
^5000 GMCs in the Galactic disk, W51 is among the top 1 % by size and the top 5-10% 
in terms of cloud mass. The 68kms~^ cloud is an elongated filament of ^ 136 pc x 22 pc 
in size. The 6:1 aspect ratio of the major and minor axis in the 68kms~^ are rare in the 
Galaxy over such a large size scale ([Solomon et al. 1987| ), and suggests that this molecular 



cloud may represent a transient feature. 

(3) The properties of the W51 and 68kms~^ clouds are compared with nearby clouds 
that have been studied in a similar manner but contain lower levels of massive star formation 
activity. While the W51 cloud is larger and more massive than nearby clouds, the mean 
H2 column density is not unusual given the large size, and the mean H2 volume density is 
comparable. The W51 GMC is similar to other clouds in that most of the molecular mass 
is contained in a diffuse molecular envelop that is not forming massive stars. The 68kms~^ 
cloud contains the largest mean column density among the clouds studied here, but this 
may be a result of an unusual viewing angle for this elongated cloud. We suggest that much 
of the star formation activity in the W51 region has not resulted from global collapse of the 
W51 cloud, but from forces acting on localized regions within the cloud. 

(4) We speculate that much of the massive star formation activity in W51 has resulted 
from a collision between the W51 and 68kms~^ molecular clouds. This conjecture can 
explain the string of embedded O stars that are spread out for 70 pc along the interface 
between the W51-68kms~^ clouds, and why massive star formation is currently confined 
to the southern ridge of the W51 GMC. 
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Fig. 1. — Galactic longitude-latitude images of ^^CO(J=1-0) integrated intensity (JT^dv) 
in 10 kms~^ intervals along the Galactic midplane < 1°) from i = 40° — 55°. The data 
are from the Massachusetts-Stony Brook Galactic Plane Survey (Sanders et al. 1986). The 
W51 GMC complex is seen most clearly in the AV = 55 — 65 kms~^ panel as bright ^^CO 
emission that extends over ^ 1° x 1° centered near £ = 49.5°, 6 = —0.2°. A high velocity 
''ridge" of ^^CO emission is evident in the AV = 65 — 75 kms~^ panel at the same £ and b. 

Fig. 2. — A 3-dimensional, velocity-Galactic longitude-Galactic latitude (^,^,6), isosurface 
plot of CO (1-0) emission at the = 4 K level for the region of the Galactic Plane 
containing W51. The data are from the Massachusetts-Stony Brook Galactic Plane Survey 
(Sanders et al. 1986). 

Fig. 3. — Channel maps of the integrated ^^CO(J=1-0) emission in 2 kms~^ intervals over 
the velocity range from 40-70 kms~^. The values printed in the upper right corner of each 
panel denote the center of the velocity intervals. The grey scales in each panel are displayed 
in a logarithmic stretch from logio(1.0 kms~^) [white] to logio(10.0 kms~^) [black]. 

Fig. 4. — Same as in Fig. 3, except for ^^CO(J=1-0). The grey scales in each panel are 
displayed in a logarithmic stretch from logio(0.6 kms"-^) [white] to logio(5.0 kms"-^) [black]. 

Fig. 5. — Top: Histograms of the number of gaussians in the ^^CO and ^^CO data cubes 
as a function of the fitted mean velocity. Bottom: The total integrated intensity in the 
gaussian components as a function of the fitted mean velocity. The designated peaks denote 
a velocity component that occurs frequently in the ^^CO and ^^CO images, and the brackets 
show the range of velocities subjectively assigned to each component. The 58 km s~^ velocity 
component is not a prominent feature as judged from this figure, but is included nonetheless 
since that is the velocity of the molecular gas associated with the brightest radio continuum 
emission. 

Fig. 6. — Various images of the W51 region. Left: Map of the ^^CO(l-O) integrated 
intensity (JT^dv) for gaussian fits with mean velocities between 56 to 71 kms~^. The 
contour levels are 10, 25, 50, 100, and 200 K kms~^. Center: Grey scale image of the 60/xm 
emission from the IRAS Sky Survey Atlas. The contours begin at logio(100 MJy/ster) with 
logarithmic intervals of 0.15. Right: Map of the A 21 cm radio continuum emission from Koo 
& Moon (1997). The contours levels are 0.015, 0.06, 0.12, 0.3, 1.0, and 2.4 Jy beam-^ The 
dashed circle represents the 0.1 primary beam attenuation level, although the image itself 
has not been corrected for primary beam response. The prominent radio continuum sources 
are labeled in the A 21 cm continuum image. The HII regions G49.4-0.3 and G49.5-0.4 are 
classically referred to as W51A, and G48.9-0.3, G49.1-0.4, and G49.2-0.4 as W51B. W51C 
has a nonthermal spectrum and is likely a supernova remnant. 
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Fig. 7. — Histogram of the ^^CO(l-O) integrated intensity for W51, the 68 km s~^ cloud, and 
clouds associated with various HII regions that have been mapped in the similar manner 
(Carpenter, Snell, & Schloerb 1995; Heyer, Carpenter, & Ladd 1996). The histograms have 
been truncated at an intensity of 3 K km s~^ since that is the highest 3a detection limit 
among the various surveys. Assuming that the ^^CO emission is optically thin, the ^^CO 
abundance relative to H2 is approximately constant, and the excitation conditions are similar, 
these histograms should reflect the relative H2 column density distributions in these clouds. 



This figure "figurel.gif" is available in "gif" format from: 
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Table 1. Properties of the W51 and 68 kms ^ Clouds ^ 



Cloud 


^min 


^ ^max 


AFfwhm^ 


Mxco Mvir 




(pc 


X pc) 


(km s-i) 


(10' Mo) 


W51 


83 


X 114 


8.9 


11.9 8.2 


68 kms-i 


22 


X 136 


4.6 


1.9 



a — Assumed distance is 7.0 kpc (Genzel et al. 1981). 

b — Full width at half maximum of the sum of the CO (1-0) gaussian fits 
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